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CT} ■ We study the interplay of neutrino oscillation and invisible decay in atmospheric 

. and long-baseline neutrinos experiments. We perform a global analysis of the full 

atmospheric data from Super-Kamiokande together with long-baseline K2K and 
MINOS in these scenarios. We find that the admixture of — > v T oscillations with 
parameters Am 2 2 = 2.6 x 10 -3 eV 2 and #23 ~ 34° plus decay of the heavy neutrino, 
z/3, with lifetime of the order r^/m^ ~ 2.6 x 10 -12 s/eV provides a reasonable fit to 
atmospheric neutrinos, although this solution becomes more disfavored (dropping to 
the 99% CL) once long-baseline data are included. Other than this local minimum, 
the analysis shows no evidence in favor of a non-vanishing neutrino decay width 
and an lower bound on the decay lifetime T3/1713 > 9.3 x 10~ n s/eV is set at 99% 
CL. In the framework of Majoron models, this constraint can be translated into 
a bound on the Majoron coupling to ^3 and an unmixed very light sterile state, 
\g s3 \ < 8.6 x 1(T 3 (2.2 eV/m 3 ). 
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Neutrino oscillations have entered an era in which the observations from under- 
ground experiments obtained with neutrino beams provided to us by Nature 
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- either from the Sun or from the interactions of cosmic rays in the upper at- 
mosphere - are confirmed and refined by experiments using terrestrial beams 
from accelerators and nuclear reactors [1]. 

In particular, with its high statistics data [2] Super-Kamiokande (SK) estab- 
lished beyond doubt that the observed deficit in the /z-like atmospheric events 
is due to — > v T oscillations, a result also supported by other atmospheric 
experiments such as MACRO [3] and Soudan-2 [4] . This was further confirmed 
in terrestrial experiments, first by the KEK to Kamioka long-baseline (LBL) 
neutrino oscillation experiment (K2K) [5], and currently by the Fermilab to 
Soudan LBL experiment, MINOS [6] 

Mass-induced neutrino oscillations are not the only possible mechanism for 
— > v T flavor transitions. They can also be generated by a variety of non- 
standard neutrino interactions or properties [1]. Prior to the highest-statistics 
SK data, some of these scenarios could provide a good description - alternative 
to Am 2 -induced oscillations - of the atmospheric neutrino phenomenology. In 
particular, it was early noticed [7] that a scenario of very fast Vj — > ui os- 
cillations plus invisible neutrino decay Vj i/jl could describe the L/E 
dependence (where L is neutrino flight length and E its energy) and the 
up-down asymmetry of the contained events in SK if sin 2 9^ ~ 0.87 and 
m j/ T i ~ 1 GeV / De (where De is the diameter of the Earth). However, with 
more precise data, it was shown that the description of the global contained 
event sample in this scenario was worse than in the case of oscillations. Fur- 
thermore for lifetimes favored by the contained event data very little conver- 
sion is expected for upgoing stopping muons in contradiction with observation. 
Based on these facts this mechanism was subsequently ruled out in its simpler 
form [8,9]. 

The possibility of atmospheric neutrino decay was revisited in Ref. [10], where 
the interplay of oscillations and decay where discussed under the assump- 
tion that oscillations where suppressed for atmospheric neutrinos, Am^ < 
10~ 4 eV 2 , so that only the mixing plus decay effects were relevant. It was 
found that a good fit to the contained and upgoing muon atmospheric data 
at that time could be obtained for Tj/m; = 63 km/GeV and sin 2 = 0.30. 
Again, this scenario became disfavored as the statistics accumulated by SK 
increased. In particular, in Ref. [11] the SK collaboration presented a study 
of the ffj, disappearance probability as a function of L/E, finding evidence 
for a dip in the L/E distribution - in agreement with the sinusoidal flavor 
transition probability predicted by mass-induced oscillations. From this they 
concluded that the mixing plus decay scenario of Ref. [10] provided a worse 
fit (by about 3.4er) than the standard oscillation hypothesis to the observed 
event distribution. 

At present, besides the strong bounds from atmospheric data, the observation 
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of the energy spectrum both at K2K and MINOS further constrains any 
flavor transition mechanism which does not lead to the correct oscillatory be- 
havior. However, this does not exclude that neutrino decay could play a role, 
even if sub-dominant, in the atmospheric and LBL neutrino phenomenology, 
and in principle affect our determination of the neutrino parameters. Con- 
versely, if this is not the case, from a joint analysis of oscillations and decay 
in atmospheric and LBL experiments one can derive a robust bound on the 
neutrino decay lifetime of the relevant states. In this paper we address these 
questions by performing a global analysis of the atmospheric and LBL data 
with — ► v T transitions driven by neutrino masses and mixing and allowing 
for neutrino decay. 

For the sake of concreteness we focus on scenarios with normal mass ordering 
of the neutrino states, m 3 > m-i > mi, in which the heaviest neutrino (z/ 3 
by convention) decays invisibly. Since v\ and V2 have large mixing with v e , 
invisible decay of these states is strongly constrained by the non-observation 
of its effects in solar neutrinos [12-14], r/m > 10~ 4 s/eV, which makes it 
completely unobservable in present atmospheric and LBL experiments. Fur- 
ther simplification arises if one assumes that the decay products are outside 
of the (u e , v^, u T ) neutrino ensemble. Indeed this is required in order for the 
decay to be fast enough because in this case the mass difference for the decay 
^3 - * ViX, Am|j, may not be directly constrained by oscillation data. This is 
the case if, for example, u 3 decays into a fourth much lighter sterile neutrino 
v s with which none of the active neutrinos mix [10]. 

With these assumptions, and neglecting the small allowed v e admixture in 
the oscillation [15], the atmospheric and LBL neutrino evolution equation 
involves only two neutrino states v 1 - = (z^, u T ). u 3 decay can be accounted for 
by introducing an imaginary part in the Hamiltonian which is proportional to 
the only relevant decay width 
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where r 3 is the z/ 3 lifetime 1 1 and U23 is the rotation matrix of mixing angle #23- 
Solving Eq. ([TJ) one gets the survival probability of v^. 

m 3 ^ m 3 L 1 /\ 777 ^ fj \ 

= cos 4 9 23 + sin 4 e 23 e'~ + 2 sin 2 6 23 cos 2 6 23 e~ ^ * cos I ^ 1 . (2) 
Equation (j2J) contains as limiting cases the scenarios explored in Refs. [7, 



1 Equation (fTJ also describes the oscillation plus decay 1^3 -> i/jl even if Vi has 
admixtures of and v T provided that its energy is degraded enough so that it does 
not contribute to the observed event rates. 
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10] (up to a relabeling u 2 <-> V3, or, equivalently, sin #23 <-> cos #23)- If one 
assumes that ArnJ, ^> E/L the oscillating term in Eq. (jSJ) averages to zero 



and = cos 4 6*23 + sin 4 #23 e T 3 B . This was the decay model proposed in 
Ref. [7]. In the alternative scenario of Ref. [10] u 3 decays into a sterile state 
Uj with which it does not mix. In this case the mass difference relevant for 
oscillations is unrelated to the mass difference between the decaying and the 
daughter neutrino states so one could have fast decays even if Am^ was very 



small. In that limit P w = y cos 2 #23 + sin 2 #23 e 2t s e J .As mentioned above 
both these limiting cases are now excluded. 

In this work we have performed a global analysis of atmospheric and long- 
baseline neutrino data in the general framework of oscillation plus decay, as 
described by Eq. (jSJ), leaving free the three parameters Am§ 2 , #23 and T3/1713. 
We have included all the SK-I and SK-II data as well as the latest K2K and 
MINOS results. Concerning the analysis of atmospheric data, an extensive 
description with all the technical details of our calculations can be found 
in the Appendix of Ref. [1]. As for MINOS, we convolve the unoscillated 
event spectrum given as a function of the true neutrino energy (which we take 
from Refs. [16, 17]) with the P MM survival probability, and with a Gaussian 
smearing function to properly account for the finite energy resolution of the 
detector. In this way we calculate the charged-current event rates, which we 
add to the neutral-current background (also taken from Ref. [17]) to obtain 
the theoretical prediction for each energy bin. As can be seen by comparing 
our event distribution for pure oscillations (lower panel of Fig. [2]) with the 
corresponding one from MINOS [17], our calculations show good agreement 
with the MINOS Monte Carlo. The theoretical predictions are then fitted 
against the experimental results, assuming a Poisson distribution with a total 
systematic uncertainty of 4%. 

Our results are summarized Fig. [TJ where we show different projections of the 
allowed three-dimensional parameter space after marginalization with respect 
to the undisplayed parameter. The hollow regions in the two lower panels show 
the allowed domains at 90%, 95%, 99% and 3cr CL from the analysis of the 
atmospheric neutrino data alone; inclusion of the LBL experiments lead to 
the full regions. The corresponding one- dimensional projections of the A% 2 
functions are shown in the three upper panels. From the figure we see that 
the best fit in this general oscillation plus decay scenario corresponds to pure 
oscillations with 



Am 2 2 = 2.4 x 10" 3 eV 2 , 9 = 45° , r 3 /m 3 > 10" 8 s/eV . (3) 



However, the figure also shows that a reasonable fit to atmospheric neutrino 
data is still possible in a oscillation plus decay scenario with 



Am 2 2 = 2.6 x 10~ 3 eV 2 , 9 = 34° , r 3 /m 3 = 2.6 x 10" 12 s/eV . (4) 
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Fig. 1. Allowed regions from the analysis of atmospheric and LBL in presence of 
— ► ^ r oscillations plus vj, decay. The three upper panels show the dependence of 
A% 2 on the parameters from the analysis of atmospheric only (dashed) and atmo- 
spheric+LBL (full line). The two lower panels show the two-dimensional projection 
of the allowed three-dimensional region after marginalization with respect to the 
undisplayed parameter. The different contours correspond to the two-dimensional 
allowed regions at 90%, 95%, 99% and 3cr CL. The lines (full regions) correspond 
to the atmospheric (atmospheric+LBL) analysis. 




This solution is well within the 90% CL regions for the analysis of atmospheric 
data alone (at A^atm = 3-8 with respect to the global best fit Eq. (j3J)), al- 
though it becomes more disfavored (at the 99% CL level, Ax ATM+LBL = 8.8) 
when LBL data are also included in the fit. The one-dimensional projection 
shows that the required lifetime at this local best fit point lies at the bound- 
ary of the 2a (3a) single parameter range allowed from atmospheric (atmo- 
spheric+LBL) data. Indeed this solution is similar to that found in Ref. [10] for 
the decay and mixing but still allowing for oscillations. Our results show that 
forcing Am^ <C 10~ 3 rules out this solution well beyond 3a (at Ax ATM = 13.7 
with respect to the global best fit Eq. (jHJ)), in agreement with the analysis of 
SK [11]. LBL data definitively rules out the mixing plus decay scenario with 
^Xatm+lbl = 39. 

A better insight on the solutions in Eqs. d3J) and (pE} can be obtained from 
Fig. [2j where we show the expected event distributions at SK and MINOS. As 
seen in the figure both solutions yield rather similar results for the atmospheric 
neutrino events. For the sake of comparison we also show the corresponding 
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Fig. 2. Zenith-angle distributions for SK /i-like events and energy spectrum at MI- 
NOS (normalized to the no-oscillation prediction) for the global best fit (full line), 
for the local best fit with oscillations plus decay (dashed line), and for a solution 
with oscillation plus decay with maximal mixing (dash-dotted line). 

distribution for a decay plus oscillation scenario with Am 2 and rjm^ as in 
Eq. (jl]) but the mixing angle still maximal. This last curve illustrates how the 
introduction of decay produces a strong deficit of atmospheric events, which 
can be compensated by the deviation of the mixing angle from maximal. Also, 
as easily seen from Eq. (f2]) this compensation is only possible with a mixing 
angle 9 < 45. The event distributions for MINOS also show how the oscillation 
plus decay scenario cannot fully account for the observed dip in the neutrino 
spectrum; this leads to the worsening of the fit when including the LBL data. 
Thus this scenario will be further tested by a more precise determination of 
the spectrum in MINOS. 



Beyond this local minimum, the analysis shows no evidence in favor of a non- 
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vanishing neutrino decay width, thus it allows to set an upper bound on the 
decay lifetime 



r 3 /m 3 > 2.9 [0.93] x 10~ 10 s/eV r 3 > 6.5 [2.0] x 10~ 10 ( m3 - \ s (5) 

V 2 . 2 (G \» / 

at the 90% [99%] CL where in the right hand side we have normalized to 
the maximum allowed value on the absolute mass scale of the neutrino from 
tritium (3 decay experiments [18,19], m 3 < 2.2 eV. 

We now turn to compare the bounds in Eq. ([2]) with the existing bounds from 
other experiments. In order to make such comparison, we must first specify 
the neutrino decay model. The reason for this is that most of the bounds in the 
literature are not derived exclusively from effects due to the neutrino decay, 
but also from effects associated with the presence of new neutrino interactions 
which are responsible for its decay. We will focus on the fast invisible Majorana 
neutrino decay u 3 — > u s J induced by the interaction 

£/ = Wa his Vj J , (6) 

where J is the Majoron (pseudoscalar) field [20-25], which has to be domi- 
nantly singlet, in order to satisfy the constraints from the invisible decay width 
of Z [26]. Alternatively for Dirac neutrinos one can have the decay channel 
v z v s x induced by a new neutrino interaction with a complex scalar field 
X [27]. In both cases the rest-frame lifetime of u 3 for m 3 ^> m s is given by 

73 = — , (7) 

where g s3 = cos # 23 g ST + sm 6*23 g sl i is the relation between the relevant coupling 
constants in the mass and flavor basis (which for Majorana neutrinos is = 
Uf a g a i3 U/3j). For these modes, the 90% [99%] bounds on Eq. (j3J) imply: 

r , 3 /2.2 eV 

g s s <4.8 8.6 xKT 3 

V m 3 

This bound can be directly compared with the constraints on the g m and g Ta 
(for any flavor or sterile state a) couplings that have been derived from their 
effect in meson and charged lepton decay. The most updated analysis [28] 
yields the model independent 90% bounds 

1^1 <9.4x 10~ 3 , |s™|<0.33. (9) 

Also, limits from decay and scattering of Majorons inside supernova yield 
bounds on \g a p\ because for large couplings the supernova energy is drained 
due to Majoron emission and no explosion occurs. However, for very large 
coupling the Majoron becomes trapped inside the supernova and no constraint 
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is possible [29,30]. As a consequence both ranges 

\g a/3 \ < 3 x 1(T 7 or \g aP \ > 2 x (10) 

are allowed [29]. 

In Ref. [31] a very strong limit was derived from the requirement that the 
neutrinos are free-streaming at the time of the photon decoupling, as de- 
duced by precise measurements of the CMB acoustic peaks, \gij\ < 0.61 x 
10 -11 (50 meV/mj) 2 . However the robustness of this conclusion has been ques- 
tioned in [32]. 

Future experiments can improve the bounds on g s/1 and g ST by orders of magni- 
tude, in particular from their cosmological effects [33] and from the observation 
of diffuse supernova neutrino background [34,35]. Till then the bounds derived 
in this work, Eqs. (jSJ) and (jSJ), from the analysis of atmospheric and LBL neu- 
trino data are the strongest applicable to the u 3 neutrino state with masses 
0(eV). In the flavor basis they represent the strongest bounds on the Majoron 
coupling to v T — u s in the full range 2.2 < < 0.05 eV allowed for the normal 
ordering of the neutrino mass states > m2 > mi with ^> m s . 
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